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Abstract 

Nonperturbative analysis of quantum entanglement and quantum teleportation pro- 
tocol using oscillator variables carried by observers in relativistic motion under the 
continuous influence of the environment is given. The full time evolution of quantum 
entanglement among static and accelerated observers is studied. The environment 
plays a dual role. While it creates bipartite and tripartite entanglement among ob- 
servers even when the initial state is separable, it suppresses the entanglement via 
decoherence. Motivated by the black hole information problem, we consider quantum 
teleportation between static and accelerated observers. Acceleration of the observer 
suppresses fidelity of teleportation. Some of the quantum information escapes outside 
of the horizon in the form of bipartite and tripartite entanglement during the telepor- 
tation process. Explicit calculation of information loss is provided. In addition to the 
loss due to the interaction with the environment, there is an intrinsic loss originated 
in a measurement process. We discuss the implications of our results on the black hole 
case. 
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1 Introduction 



The equivalence principle tells us that a particle under constant acceleration can be viewed 
as a particle under the influence of a static gravitational fleld. Explicit coordinate transfor- 
mation shows that a particle moving in a hyperbolic trajectory due to constant acceleration 
in Minkowski spacetimc is equivalent to a particle moving in a geodesic in Rindler spacetime. 
Under the coordinate tranformation, positive frequency modes in Minkowski spacetime are 
expressed as the mixture of positive and negative frequency Rindler modes. This implies 
that the Minkowski vacuum annihilated by positive Minkowski modes is not equivalent to 
the Rindler vacuum. As shown in [1], an accelerated observer sees the Minkowski vacuum 
filled with Rindler particles. The spectrum of Rindler particles is thermal, similar to the 
spectrum of the Hawking radiation from a blackhole[2]. The Unruh effect modifies the 
quantum dynamics of accelerated systems in a nontrivial way causing apparent discrepancy 
between the interpretation of radiation emission and absorption between accelerated and 
static observers[3]. 

Stimulated by recent progress in quantum information science, various attempts to un- 
derstand relativistic quantum phenomena from the view of quantum information have been 
madc[4]. Quantum tclcportation between a static and an accelerated observer was discussed 
in [5], where one of the cavities containing a single mode is assumed to follow the accelerated 
trajectory. Because there was no consideration of dynamical aspects and other crucial fea- 
tures in [5], their derivation is not expected to be applicable in the general situation[6]. The 
entanglement in the system with one accelerated and one static[7] and both accelerated [8, 9] 
was discussed. Entanglement dynamics involving accelerated observers was studied in [10]. 

On the problem of black hole information loss[ll], entanglement in the Hawking pair is 
often regarded as a crucial ingredient to rescue information from a black hole in the form 
of the Hawking radiation [12]. The nonlocal process present in quantum teleportation is also 
considered to play a role[13]. More recently, it was speculated that if the collapsing matter 
and incoming Hawking particle at the black hole singularity [14] are maximally entangled, 
information can escape from a black hole as a process similar to the quantum teleportation. 
This scenario, unfortunately, is likely fiawed by the interaction between the collapsing matter 
and the incoming Hawking fiux[15]. In the case of the qubit teleportation assuming the 
random unitary interaction between the original state and one of the entangled pair [16], 
however, nearly perfect information can be restored with proper encoding after quantum 
error correction. It is of interest to see if the same scenario works in the spacetime with 
a horizon. Motivated by these works, in this paper, we study the dynamical evolution of 
quantum entanglement and fidelity during quantum teleportation due to accelerated motion 
of observers. Near horizon, the Rindler metric is equivalent to the Schwartzchild metric. 
Thus our work can be viewed as performing a quantum simulation of information fiow around 
the black hole horizon. 

First we develop a general framework to study the entanglement among an arbitrary 
number of observers generalizing the formahsm in [25] to the relativistic setting. Each 
observer's quantum state is assumed to be described by a harmonic oscillator, which couples 
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linearly with a scalar quantum field. In this work, we consider the vacuum state of the 
field as environment. Wc are primarily interested in the quantum state of observers and 
trace out all the field modes in our calculation. Then we study the entanglement dynamics 
among observers' quantum states when each observer is allowed to move in a prescribed 
trajectory. When one of the observers is subjected to constant acceleration, his worldhne 
in a Minkowski coordinate becomes a hyperbolic curve and asymptotically approaches his 
event horizon. The observer sees the Minkowski vacuum as a thermal state. Entanglement 
between his and other observers' states will be modified accordingly. 

Perturbative approximations such as Born or Born-Markov approximations are widely 
used for the study of decay processes and estimation of various quantities of open systems. 
Radiation processes from the accelerated particles and near black holes [17] are often stud- 
ied under these approximations. However, naive perturbation expansions do not guarantee 
the positivity of the reduced density matrix. Entanglement properties are measured by the 
negativity of the partially transposed density matrix and sensitive to the positivity of the 
density matrix. For this reason, in this paper, we try to use an exact expression without 
making perturbative approximations. Nonperturbative analysis is one of the key achieve- 
ments in solving the apparent paradox in the radiation-reaction process viewed in inertial 
and accelerated frames [18]. 

Fluctuations of the canonical phase variables obey uncertainty relations. Uncertainties 
are modified in the presence of the environment. In the system out of equilibrium, they 
are time-dependent. In order to set up the quantum teleportation protocol [19] using our 
oscillator variables, wc introduce another observer who carries a state to be teleported. 
We study uncertainties in the niTilti-modes and follow their dynamical evolution. We are 
particularly interested in the case when one observer is under constant acceleration and study 
how motion of the observer affects uncertainties. Although uncertainty relations should 
always be satisfied by any physical states in principle, they are no longer guaranteed under 
perturbative approximations. Similarly to entanglement, they are sensitive to the positivity 
of the density matrix. Next we study entanglement among these observers. As shown 
in [20], uncertainty relations in the partially transposed density matrix give the necessary 
and sufficient condition for separability. Entanglement measures, negativity and the log- 
negativity, will also be calculated from the uncertainties. 

Next we study entanglement among these three observers and clarify the relation between 
bipartite and tripartite entanglement. For a three- mode oscillator system, it is possible to 
find the genuine tripartite entanglement Eabc from the monogamy relation[21, 22]. In the 
conventional quantum optical setup for the teleportation[23], the effect of vacuum is treated 
locally and there is no role for the tripartite entanglement. In our case, vacuum is a ground 
state of the quantum field which spreads nonlocally in spacetime. Therefore its infiuence on 
the observers will develop nonlocal correlation between their states. In order to see this, we 
set the initial state in a product state between a coherent and a two-mode squeezed state 
and see how the system develops tripartite entanglement dynamically in time. 

Then we perform quantum teleportation using an entangled pair one of which is carried 
by a receiver under constant acceleration. We include the effect of quantum field continuously 
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interacting with the system and calculate time-dependent teleportation fidelity as a measure 
of success of teleportation. First we study the case with arbitrary final states of the sender 
and the carrier of one of the entangled pair by summing over all possible states between them. 
Next we consider the generalized final state to which their state will be projected after the 
measurement to be in a two-mode squeezed state. We send the measurement result to the 
receiver followed by the appropriate unitary transformation on his state so as to recover the 
original state with generic initial and final squeezed states and see how niTich information 
can be recovered. We consider both in the presence and absence of the interaction between 
the system and field. 

We first obtain the general expression for the fidelity for arbitrary amount of squeezing for 
the initial and final state. Then we look at the limiting case when the final state is maximally 
entangled. In the conventional quantum teleportation scheme, the final state is generally 
assumed to be maximally entangled while the initial squeezing is finite due to experimental 
limitation. Since our motivation is to probe black hole information loss from the possiblly 
unknown final state, we consider the generalized final state by introducing a finite squeezing 
parameter. In all cases, we obtain the exact dynamical expressions for the fidelity. We also 
study the asymptotic long time limit of the fidelity and discuss the implications of our results 
on the black hole information problem. 

In Sec. 2, we develop general formulation for N observers each carrying quantum state 
moving in a four-dimensional spacetime interacting with quantum fields. Tracing out quan- 
tum fields, we obtain the reduced density matrix for H observers' quantum state. In Sec. 
3, from a correlation matrix of the system, we calculate uncertainty relations for multi- 
modes. Also from the correlation matrix after partial transpose, we calculate the bipartite 
entanglement measures for two and three modes from which we obtain genuine tripartite 
entanglement. We give a detailed study of the time evolution of these quantities under the 
infiuence of the environment in the absence and presence of acceleration of one observer. In 
Sec. 4, quantum teleportation among three parties is studied. We first study the case by 
integrating out all the final states and look at the time evolution of the teleportation fidelity. 
While introducing the basic scheme of teleportation, we generalize the possible final state to 
an arbitrary two-mode squeezed state. Then we consider the effect of the environment for 
the protocol and study the time evolution of the fidelity. We obtain the result for the max- 
imally entangled final state by taking the infinite squeezing limit of the final state. Finally 
we consider the maximally entangled initial and final states including the environment. We 
will study the full time evolution and the long time, weak coupling limit of the fidelity. We 
will use natural units: G = h = c = kB = l- 

2 General formulation 

We consider the system composed of observers moving in prescribed trajectories. Our 
Hamiltonian is a relativistic generalization of a jV-Brownian oscillator model studied in 
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[24, 25] and given by 

Hs^J:^ + Vo(R,,...,Rm), 

3=1 ^ 

where the bare potential Vq is a sum of the physical potential 
and the counter term AT/ 



N N AT/ 

j=i 1=1 ^ 



where 



TT 

with the frequency cutoff A for the scalar field below. 

Each oscillator variable Rj couples linearly with the scalar field 0j at Xj as 

M 

Hi = I^Aj- i^j(r,)0j(xj(r,),i(rj)). 

3=1 



We assume that Rj is an internal coordinate of a j-th observer which depends on a 
time Tj of the observer. Our total Hamiltonian is thus Lorentz invariant. 
The field Hamiltonian for each (f)j in a D-dimensional space is 

1 

2 



3^3 



Each field 0,- allows a mode decomposition: 



with Uk = vW+m^. 

The Heisenberg equations that Rj satisfy are 



MjRjir) + Mjn]Rj + ^VjiRj + 2 E dsaiji(Tj, si)Ri(si) = 0, 

I I ■'^ 



where 



(^ijiiTj, t[) = -X] COS [k{xj{Tj) - Xi(r/))] sin [uJk{t{Tj) - t(r/))] /2uJk (9) 

k 

is an imaginary part of the response function [28] defined as 

a,,(T,-, r/) = \] Ek cos [k{x,{T,) - xz(r/))] e'^^^^^^^^ /2^fc, where Ek = I rf^A;/(27r)^. We 
set the initial condition at t = hypersurface when all observers are at rest with their proper 
times set to tj = for all j. Eq.(8) contains nonlocal kernels aiji{Tj,T[) and the exact time 
evolution of Rj at r depends on its and other observer's past history. Thus the dynamics of 
i?i is not Markovian. 




Figure 1: The world line of three observers. 

Let us consider the case A/" = 3 in which the first observer Bob is accelerated with 
a constant acceleration a in an accelerated trajectory so that its coordinate is given by 
t(r) = sinh(ar)/a and Xi{t) = cosh(ar)/a. The second observer Alice and the third observer 
Chris's locations are fixed at X2 = and 0:3 = —L < 0, respectively. We assume that Alice 
and Bob are close to each other (e << 1) and share a common environment (02 = ^s), 
while the accelerated observer Bob is located in a distance such that -Ri can be treated as 
independently coupling to an environment. The response functions (9) are then given by 

an(r,r') = A? cos [A;(xi(r) - xi(r'))] e-*-^(*M-*(-')V2a;fc, 

k 

a,,{r,r') = ^ e-^-'=(^-^'V2^fc (for j = 2, 3), 

k 

a23(r,r') = a32(r, r') = A2A3 E cos(Lfc)e-*'^'=(^-^')/2^fc, (10) 

k 
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and ai2 — an — 0(21 — Oi^i — 0. Note that the coordinate time t and the proper time are 
the same for Ahce and Chris but not for Bob. The most interesting term above is q;ii(t) 
where each mode in Minkowski spacetime is expressed in Rindler coordinates. The conver- 
sion of modes in different spacetimes is given by the Bogohubov transformation between 
annihilation/creation operators [27]. Writing the momentum space integral in (7) in terms of 
Rindler modes reproduces a standard derivation of the Unruh effect based on the Bogohubov 
transformation[27]. For a massless field 0i^2 in ^ three dimensional space, the kernels in (10) 
can be brought to the following form: 



cos(A;t) coth(— ) — isin(A;r) 
a 



"^^^^^ = 2^-1 

/-A 

'^jji'^) ~ J dkk [cos{kT) — i sm(kT)] , (forj = 2,3), 

A2A3 /"A 

a23(T) — Q;32(t) = — — - / (i/c/csin(LA;) [cos(A;t) — isin(A;T)] . (11) 

vr^L Jo 

The kernel 0:23 is responsible for the interaction between Alice and Chris's state that becomes 
negligible for large separation as 1/L, which is a characteristic decay of a massless field in 
three dimension and justifies neglecting the interaction between the accelerated observer and 
the rest of observers. Comparing the first kernel an with that of the infiuence functional[28], 
we see that the first observer sees the vacuum as a thermal bath with the Unruh temperature 
Tu = a/27r[29, 30]. Writing the influence kernel directly in the Rindler time has advantages 
in studying quantum dynamical aspects of observer's quantum states under the influence of 
the Unruh effect, which we will describe below in details. 

For large A and small separation L, the Heisenberg equations of motion for three coor- 
dinate variables can be written in the local form: 

Miiiiit) + MinlRi{t) + r,Ri{t) = 0, (12) 

M2i?2(t) + M2nlR2{t) + 724(t) + liisit) = 0, (13) 
MsRsit) + M3QlRsit)+^sMt)+lMt) = 0, 

where we write 7^- = A^/tt (for j = 1,2,3) and 7 = X2X3/tt. The divergence from nonlocal 
terms are canceled by counter terms. For our case, AV12 = AV21 = AVi^ = AV31 = 0. 
Hereafter we will set Mi = M2 = M3 = 1. 

We write the solution for the first equation in (12) with initial conditions i?i(0) = and 
i?i(0) = las5i(T): 

9.ir) - ^^e--, (14) 



where = — 7^ . A pair of solutions of (13) with initial conditions -R2(0) = -R3(0) = 
and i?2(0) = 1, -^3(0) = will be written as h-i{T) and h^ir), respectively. For Q2 = ^3 and 
A2 = A3, the solutions are given by /i3(t) = (5(3 (t) + go{t))/2 and /i5(t) = (5-3 (t) - go{t))/2, 
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where 



2r 



(15) 



are the solutions corresponding to two normal modes of a coupled oscillator and fl^^ = 



General solutions with arbitrary initial conditions Rjo and Pjo of the Heisenberg equations 
(12) and (13) for j = 1,2,3 are 



CRjR^RkO + CRjPkPkO + / dsg2i^-i{T - s)0j(s), 

k=l k=l ''^ 

3 3 



Pji^) = Y Cp.R^Rko + Y CpjPkPko + Aj- / dsg2i^ (r - s)0j(s), 
k=i k=i ''^ 

where g2j = g2j-i and h = 1, l2,3 = 2. 

The expectation value of phase space variables can be expressed in the matrix form: 



(16) 







/ -Rio 


\ 


1 Cr^r-^ 


CriPi 














\ 


( Rio 




(Pi) 




Pio 






CpiPi 
















Pio 




{R2) 




R20 










CR2R2 


CR2P2 


CR2R3 


CR2P3 




R20 




{P2) 


P20 










CP2R2 


CP2P2 


CP2R3 


CP2P3 




P20 








R30 










CR3R2 




CrsRs 


CR3P3 




R30 




I (Ps) ) 




\ P30 


J 


V 





CP3R2 




CP3R3 


CpaPa 


J 


\ P30 


J 



(17) 



The time evolution matrix C for our solutions is given by 



f fi gi \ 

/2 ^2 

0/3/^3 h h 

0/4/14 /e h 

0/5/15 fs h 

V h h h J 



C = 



;i8) 



where /aj = hj-i, and /a^-i = g2j - 2g2j{0)g2ij-i for j = 1 and /a^-i = /i2j - 2h2j{0)g2ij-i 
for j = 2, 3. 

It is convenient for our purpose to use a characteristic function for the Wigner representation[31] . 
For the AA-particle system: 



Xw{y,r) = IV 



p(0)e 



(19) 



where we defined X2j_i = ^JVljRj, X2j = Pj / ^Vlj ioi j = 1, 2, A/'and 3^ = {Yi...Y2jKf)- The 
matrix C is scaled accordingly. Now we trace out the field in order to obtain the reduced 
dynamics of the system. In our case J\f = 3 with a factorized initial condition: p(0) = ps(0) ® 



8 



nLiP^,(0), xw is also factorized into two parts as Xw{y,T) = Xw'l^^ IlLi xfe(3^, t). 



The system part Xw{y^^) = Tr^ 



P5(0)e 



, where Xcjir) are solutions of the 



Heisenberg equations with (pj — 0, and the environment parts 
Xw(y^^) = Tr^i P<^i(0)exp 
xtviy^^) = Tr^2 P<^2(0)exp 



(20) 



where [Z] is equal to 1/2 (for even I) and (Z + l)/2 (for odd I) and (Z) = Z for Z = 3, 4 and 
(Z) = / — 2 for / = 5, 6. In response to the state of the system, each field mode is shifted 
by system-bath interaction. After the time scale of each field mode, the back action of this 
process changes the system's state depending on the shift of the field mode. The collective 
effect of these processes leads to the non-Markovian evolution of the reduced system. 

The field characteristic function in (20) can be evaluated exactly. We assume the envi- 
ronment is initially in a Minkowski vacuum state. Its density matrix is given as ^^(0) =| 
Om) (Om I ■ The environment characteristic function can be written as 



xt^{y,t)xUy,t) - exp 



exp 



--y^^y 



(21) 





S12 














\ 


S21 


S22 






















S33 


S34 


S35 


^36 












E44 


S45 


^46 
















^56 




I 





^63 


^64 


^65 


^66 


/ 




where 



S,7(r) = 



3 S lo Jo ds'gjir - s)gi{T - s')e''=(^iW-^i(^'))e-i'^fc(*i(^)-*i('*')) for j, / = 1, 2 



are time-dependent, nonequilibrium fluctuations of the system variables induced from the 
environment. Nonvanishing off-diagonal correlations Sjj for j ^ I generated from the inter- 
action with the environment are signatures of induced correlation and entanglement among 
observer's quantum states. For j, / = 1,2, 

^ni^) — ~ dujujcoth( — ) ds ds'gjlr — s)gi(T — s') cosuls — s'). (22) 
2tt Jo a Jo Jo 

If the initial states of system variables are all Gaussian states with vanishing mean posi- 
tions and momenta, {X{0)) — 0, the system characteristic function also takes the Gaussian 
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form: 



exp 



(23) 



[ {{Xce,Xci}) 




where {A, B} = (AB + BA)/2 is an anticommutator and Xc = {Xc'i---Xce) satisfy the 
equations of motion (12) and (13) for damped harmonic oscillators. (AA')^(r) are essentially 
the initial fluctuations of the system variables shifted by the damped oscillatory motion of 
a coupled harmonic oscillator. Combining with the characteristic function for the field, we 
obtain 



Xw{y,t) = exp 



exp 



= exp 



{{X,,X,}) ... {{X,,X,}) 



(24) 



{{Xe,X,}) ... {{X,,Xe}) J \Ye 
/ {{X,c,X^c}) + ^u ... {{Xci,Xc6}) + ^i6\ f Y^\ 



\ {{Xc6, Xci}) + '^61 ... {{Xc6, Xce}) + ^66 J 



\YJI 



3 Entanglement dynamics of accelerated oscillators 

A separability criterion for a bipartite two-level-system can be naturally extended to contin- 
uous Gaussian variables [20]. The necessary and sufficient condition for separability of the 
density matrix is to have only non-negative eigenvalues after the partial transpose of one 
of its subsystem. For Gaussian variables in the Wigner distribution, the partial transpose 
of a density matrix in one of the oscillator component is equivalent to a mirror reflection 
of that component. The necessary and sufficient condition for the general (1 + N") bipar- 
tite Gaussian modes to be separable is that the partially mirror reflected state is still a 
physical quantum state that satisfies the uncertainty principle [32]. For instance, for three 
oscillators, the bipartite entanglement between the first and the rest of oscillators can be 
measured by first taking a partial mirror reflection on the first variable in the phase space as 
(i?i. Pi, i?2, -P2, -R3, -P3) — {Ri, —Pi: R2, P2, -R3, -Pa), then seeing if the resulting state satisfies 
or violates the uncertainty relation. Note that the total mirror reflection does not change 
the entanglement properties thus the reflection on the second and third variables yields the 
same result. In terms of X, this can be expressed as a matrix multiplication by the matrix 
77 = diag(l, —1, 1, 1, 1, 1) as X ^ r]X. It follows that the partial mirror reflection transforms 
the covariance matrix as 



2„r 



{i^xy ^ 7]{AXY7] 



(25) 
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Prom Williamson's tlieorem[33], there exists a symplectic transformation that diagonal- 
izes any positive-definite 2n x 2n symmetric matrix into the following form: 













\ 





Ci 














Cn 







I 







Cn 


I 



(26) 



Although such a symplectic transformation docs not preserve the eigenvalue spectrum in 
general, the diagonal components Q for / = l.-.A/" can be calculated as follows. Writing a 
commutation relation in a 2J\f x 2J\f matrix form as [X,, Xj] 



iVij with 



/ 1 
-1 



\ 




V 








1 





we construct a real symmetric matrix /S.XT [AXYV^ AX . This matrix has an eigenvalue 
spectrum C,f {I = l...A/')[34]. The uncertainty relation can be generalized to a symplectic 
invariant form {AX^ + ir/2 > 0. By changing to the diagonalized form (AA"^)^ above, the 
uncertainty relation is equivalent to saying that > 1/2 for all /. 



3.1 Bipartite entanglement 



For a bipartite continuous variable system, there is a criteria for a necessary and sufficient 
condition for separability [20, 35, 36]. We write the 2x2 covariance matrix here as 



(27) 



then eigenvalues (f can be written explicitly in terms of the following symplectic invariants 
Ai_2 constructed from the determinants of covariances 1^41, |£'i,2|, K^'^)^! as 



Then 



Ai 

A2 



|Di| + ID2I +2|A|, 
\{AXn 



Ai ± JAf - 4A2 



Under the partial transpose (25) without the third varible, A 
will be changed to 



Ai±7Af-4A2 



(28) 

-A and these eigenvalues 
(29) 
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where Ai = \Di\ + \D2\ — 2|A|. The separabihty conditions are 

XI > i (30) 

Note that the inequahties for and A_|_ follow automatically from those for C_ and A_. 
Thus A_ carries the essential information on the separability of quantum states. 

Let us consider a two-mode squeezed state with a squeezing parameter r as an initial 
state[37]. Its correlation matrix is 

In the Wigner representation, the same state can be expressed as 

W(Ri, i?2, Pi, P2) = Ae-e"^K«i-R2)2+(Pi+P2)Vf^]-e-2'^K«i+fi2)2+(Pi-P2)Vf^]. (31) 

This state can be obtained by acting a squeezing operator e*'"(^i^2-PiR2) ^y^q vacuum. For 
a large squeezing r — 00, W{Ri, R2, Pi, P2) ~ 5{Ri — R2)5{Pi + P2), thus it becomes the 
EPR (Einstein- Podolsky- Rosen) state [38]. 

In Fig. 2, the temporal behavior of the uncertainty C- as a function of the proper 
time of Bob is plotted. The initial state is a pure two mode squeezed state introduced 
above. This state satisfies the uncertainty relation with the minimum uncertainty 1/4. As 
the state becomes mixed, the uncertainty first increases in time and oscillates periodically. 
Acceleration of the observer yields larger uncertainty by making the amplitude of oscillations 
larger. In Fig. 3, the temporal behavior of C- (uncertainty) and A_ (uncertainty after the 
partial transpose) is plotted. The initial state is a pure two mode squeezed state with r = 0.5. 
The state exhibits environment- induced disentanglement at Vtt ~ 0.4 with no acceleration 
and at Vtt ~ 0.31 with acceleration. The role of acceleration is similar to that of temperature 
that increases the disentanglement rate. 

Now we study entanglement measures; the negativity and the log negativity. The nega- 
tivity j\/'[39, 40] for our system can be defined as 

Wp^W - 1 

U = (32) 

where is the reduced density matrix after the partial transpose. M is equal to the sum of 
all negative eigenvalues of . Thus it measures how much fails to be positive. The Peres 
criteria[35] tell us that this can be used as a measure of entanglement. It is an entanglement 
monotone, which does not increase under local operations and classical communications. 
The logarithmic negativity Ej^ defined as 

£V = ln||p^|| (33) 

is also an entanglement monotone. The diagonalization of (AA")^ changes the original state 
into the thermal state. Then the partially transposed density matrix p^ after the same 
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Figure 2: Uncertainties ( C- in Eq. (28) ) in the presence and absence of Bob's acceleration 
are plotted. The initial condition is a two-mode squeezed state with r — 0.1. fl is set to 1 
for all numerical plots. Other parameters are 7 = 0.001, A = 100. 



transformation also has the thermal form. It can be written in terms of the symplectic 
invariants X± as 



Pr 



- n 




n±){n±\ 



(34) 



For separable states, X± > 1/2. Then ||p^|| = 1 and jV = Ej^ — 0. For entangled states, 
A_ < 1/2 but A+ > 1/2. Thus both J\f and E/j- can be expressed in terms of A_ as 



Em 



max 



0, 



1-2A_ 



4A_ 

max [0, - ln(2A-)] . 



(35) 



In Fig. 4, the negativity A/" and the logarithmic negativity Ej^ are shown as a function of the 
proper time. The initial state is a two mode squeezed state. They both vanish at the same 
time since they both give the necessary and sufficient condition of bipartite entanglement. 
The larger acceleration yields the larger disentanglement rate. 



13 




Figure 3: The temporal evolution of uncertainties ( C- in Eq. (28) and A_ in Eq. (29) ) 
before and after the partial transpose in the presence and absence of Bob's acceleration is 
plotted. The initial condition is a two mode squeezed state with r = 0.5. Other parameters 
are 7 = 0.05, A = 50. 



3.2 Tripartite entanglement dynamics 

It is convenient to characterize an n-mode uncertainty relation as a sum of symplectic 
invariants[32] as 

E„ = X:4^-"(-l)"+^AJ, (36) 
j=o 

where A!?(j = 1, n) are principal minors [41] of the matrix r(AA')^ of order 2j. We defined 
A^ = l. 

For the matrix in the Williamson normal form (26), A" = J2{ki,...,kj} Cfci---Clj(j = 1) ^)) 
where {ki, kj} is an unordered set of different integers between 1 and n. Then we have 

= 4^(1- E < + ••• + (-4r E C.V-CL) (37) 

fcl = l fel = ...=fen=l 

Thus the uncertainty relation implies E„ > 0. 
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Figure 4: The temporal evolution of the negativity (32) and the lognegativity (33) for dif- 
ferent acceleration before and after the partial transpose is plotted. The initial condition is 
a two mode squeezed state with r = 0.1. Other parameters are 7 = 0.1, A = 50. 

For the tripartite case n = 3, we write 



(38) 



Then 



A? 



\Di\ + ID2I + l^^al + 2IA12I + 2|Ai3| + 21^231 

+ |L>2| + IDs] + 2|ii2| + 2|ii3| + 2|i23| 

l(AA')^ 



(39) 



where "stands for the codeterminant. In Fig. 5, E3 is plotted as a function of time. E3 > 
throughout the entire evolution indicates that the uncertainty relation is always satisfied. 
The effect of acceleration is again similar to the effect of temperature which increases the 
uncertainty. 

Under the partial transpose on the first variable, |yli2|,|Ai3|,|/li2|,|/li3| change signs and 
the uncertainty after the partial transpose E3 becomes 



1 A? 

E3 = + ^ 

^ 64 16 



A3 
^2 



+ A^, 



(40) 
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Figure 5: The temporal evolution of the three-body uncertainty S3 is shown. The initial 
condition is a two mode squeezed state with r = 0.2 and a ground state of Chris. Other 
parameters are 7 = 0.005, A = 50. S3 > indicates that the uncertainty relation is always 
satisfied. 



where 



A? : 
A3 : 

4 = KAA") 



\Di\ 


+ 




+ 


\D3\ 


-2IA12I 


-2|Ai3| + 2|A23| 


\Di\ 


+ 




+ 




-2|ii2| 


-2|ii3| + 2|i23| 



Prom (37), 



S3 = {xi-\)-{>^l-\)- (41) 



For separable states, each term on the right hand side is nonnegative, thus S3 > follows. 
For our interest of (1 + 2) mode bipartite system, or more generally, (1 + n) bipartite system, 
> is shown to give a necessary and sufficient condition for separability [32]. 

Three-mode entanglement for continuous variables can be categorized into the following 
sets[42] based on the values of the bipartite entanglement measure Ei with i = A{BC), {AB)C, {AC)B 
{Ea(bc) is bipartite entanglement between Ahce and Bob-Chris together): (1) fully entan- 
gled states (none of is vanishing.), (2) one-mode biscparablc states (only one of Ei is 
vanishing.), (3) two- mode biseparable states (only two of Ei arc vanishing.), (4) three-mode 
biseparable states (all of Ei are vanishing but these states cannot be written as a mixture of 
tripartite product states), (5) fully separable states (mixture of tripartite product states). 

For bipartite entanglement for three qubits, the square of the concurrence C is known to 
satisfy the following monogamy inequality [21] 

^A{BC) — ^AB + ^AC^ (42) 
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where Cab is the concurrence between A and B, etc. The inequahty also holds for Gaussian 
states if we take the square of the logarithmic negativity [22] . Genuine tripartite entanglement 
for continuous variables Eabc can be defined by minimizing the difference between the left- 
hand and the right-hand side for all possible combinations of three subsystems, namely, 

Eabc = minj^-fe - E^j - E^jfej . (43) 

Note that the monogamy inequality guarantees Eabc > 0. In Fig. 6, the temporal evolu- 
tion of bipartite entanglement Eab,Eac,Ebc arc shown. The initial state is chosen to be 
the product state of a two-mode squeezed state of Alice and Bob and a coherent state of 
Chris. Thus initially Eac = Ebc = 0. As time evolves, decoherence due to the interaction 
with the environment [43] causes Eab to damp out. The interaction with the environment 
induces entanglement between Alice and Ghris's state, which evolves periodically in time and 
some of entanglement in Eab tranfers to Eac- There is no entanglement between Bob and 
Chris's state during the entire period since there is neither initial nor interaction-induced 
entanglement between them. Note a dual role of the environment. It is both the source 
of multipartite entanglement between observers and of decoherence that decreases entangle- 
ment. The former causes initial increase and the latter causes damping of entanglement in 
time. The similar behavior is already seen in two static Brownian oscillators[25]. In Fig. 
7, the temporal evolution of Eabc is shown. The initial state is the same as in Fig. 6, a 
product state of a two-mode squeezed state of Alice and Bob and a ground state of Chris. 
Thus initially there is no genuine tripartite entanglement. The nonlocal vacuum generates 
entanglement among three variables. The same vacuum also induces decoherence which 
eventually washes away all entanglement. The Unruh effect suppresses both bipartite and 
tripartite entanglement. 




Figure 6: The temporal evolution of bipartite entanglement Eab, Eac, and Ebc are plotted 
as a function of time. The initial condition is a product state of a two-mode squeezed state 
of Alice and Bob with r = 0.3 and a ground state of Chris, a = in Fig.6(a) and a — SOtt 
in Fig. 6(b). Other parameters are 7 = 0.001, A = 50. 



17 




' 1 ' 1 ' 1 ' 1 ' 1 ' J — I — ■ — I — ■ — I — ' — I — ■ — I — ' — 1- 

0.0 0.2 0.4 0.6 0.8 1.0 2 4 6 8 10 



Figure 7: The temporal evolution of the genuine tripartite entanglement Eabc is plotted as 
a function of time. The initial condition and unshown parameters are the same as in Fig. 
6. The short time behavior is shown in (a). While squeezing of the initial state enhances 
Eabc-i acceleration suppresses it. The long time behavior in (b) is damped oscillatory. 



4 Quantum teleportation with an accelerated observer 

Many interacting quantum systems can be simulated with a proper design of quantum 
circuits[44], where the interaction among subsystems are mimicked by quantum gates. Along 
with the properly chosen measurement and post-selection scheme, we can simulate many in- 
teresting systems under nontrivial situations. In [14], similarity between information escaping 
from a black hole and a quantum teleportation protocol was pointed out. In the presence of 
the generic interaction between collapsing matter and a Hawking particle, however, unitarity 
cannot be completely rescued[15]. Encoding quantum states with quantum error correction 
may still recover the full information. 

In these works, only speculations about the possible quantum process in and around 
a black hole are given based on the analogies with quantum circuits. No analysis on the 
curved spacetime was provided. It is important to note that qubits studied there are in a 
flat spacetime without event horizon while information loss is predicted in quantum theory 
in curved spacetimes with event horizon that hides information. 

Since our model can simulate the quantum dynamics of many observers moving in an 
arbitrary fashion, it provides a good setting to study the Unruh effect on various quantum 
protocols. From the equivalence principle, the acceleration due to gravitational force can be 
locally transformed away by the appropriate coordinate transformation. In the absence of 
other forces, the observer then follows a geodesic in a curved spacetime. The Unruh effect 
on a constantly accelerated observer can be identified as the effect on the static observer 
in Rindler spacetime. Since near the Rindler horizon, the Rindler spacetime looks similar 
to the Schwartzchild black hole spacetime, this is equivalent to studying the Hawking effect 
on the quantum dynamics of a static observer near a black hole event horizon. In Fig. 8, 
our quantum teleportation scheme in (b) is compared to the scenario in [14] of information 
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Figure 8: Penrose diagrams of (a) Schwartzchild spacetime with a falling matter and a 
Hawking pair and (b) Rindler spacetime with our quantum teleportation scheme. Shaded 
region contains a pair of wormholes through which we can retrieve measurement results 
beyond the Rindler horizon. 
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escaping from a black hole in (a). In both figures, two parties are falling inside of the horizon 
and hit one region while the other party is moving toward future without falling inside of 
the horizon. The role of singularity in (a) is viewed as an effectively boundary condition, 
while it is replaced by the measurement performed at some time in (b). If Bob is under 
constant acceleration indefinitely, once Ahce passes the event horizon, we need an additional 
mechanism to send the measurement result to Bob for the teleportation to be successful. In 
the figure, a pair of wormholcs arc used to retrieve the measurement result. In [14], this 
was simply assumed as pointed out in [15]. In the present work, we view our scheme as a 
simulation of the process possibly occurring around a black hole and are only interested in 
the process until the measurement is performed. Bob can stop accelerating eventually and 
receive the message from Alice. This is indeed the situation analogous to a decaying black 
hole[45]. Although our system docs not include gravitational interaction between observers, 
as we saw in Sec. 3, the background field induces the effective interaction and entanglement 
between them. 

Quantum teleportation with continuous variables based on [23] have been demonstrated 
experimentally [46]. Initially Alice and Bob share a two-mode squeezed state. The mode 
carried by Alice is mixed with Chris's state prepared in a coherent state through the 50 
% beam-splitter. Their state is measured by the homodyne detector and the result is sent 
to Bob through the classical channel. Bob's state is shifted according to the measurement 
result in order to reproduce the same state as Chris's. 

Now let us first discuss our problem, instead of measuring the state of Alice and Chris, 
by tracing over all states for Alice and Chris at some time in the future[28]. Since we are 
concerned with the Unruh effect on the quantum teleportation protocol, first we consider 
the influence of the environment only on Alice and Bob's state. Then Chris's state remains 
in a pure state and the density matrix is the direct product of Chris's and Ahce and Bob's 
state as \oi)cPAB{o(\c, where pab can be a mixed state density matrix. The fidelity F is an 
overlap between Bob's final state and Chris's initial state given by 

(44) 



where Vf^ = (AA')?, - {AX)l2 - 2iAX^^ and Vf^ = {AX)j^ + {AX)l^ + 1 and = 
■p(o)2 _ |p(o)|2 j|- jg gg^gy ggg ^Yiai the fidelity for this process gives an one-mode Q 
distribution function for Bob's state. The more Alice and Bob's quantum state is squeezed, 
the more it is entangled. Then after tracing out Alice's degree of freedom, Bob's state is more 
mixed and the less fidelity results. This is opposite to what we expect in the teleportation 
scheme where large entanglement between Ahce and Bob yields high fidelity in the teleported 
state as we will see below. 

The fidelity takes its largest value at the center a = as F„iax = 20^^^^^^"^ . After a long 
time, Bob's state thermalizes at the Unruh temperature. The off-diagonal components of 
the correlation matrix will vanish and, in the weak coupling limit, their diagonal components 
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F = 



B{a\TrAc 
2 



W)cPab c(tt|] 



2lf) 
L>(o) 



are given by {AX)l^ = {AX)l^ = coth(7rQ/a)/2. This gives D^°^ = (coth(7rQ/a) + 1)^ and 
the long time hmit of F can be written as a function of acceleration as 



F(oo) 



■ exp 



coth(7rQ/a) + 1 



\a\ 



(45) 



coth(7rQ/a) + 1 

We see that this is independent of the initial squeezing. The peak value F^axi^) — 
2/ [coth(7rf2/a) + 1] takes from 1 to as the acceleration a varies from to cxd. We see 
that Fmax measures the overlap between Bob's and the ground state. For the zero acceler- 
ation case, a = 0, Fmax{oo) = 1, thus there is no loss of information. In Fig. 9, the time 



(a) 
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Figure 9: The temporal evolution of the fidehty in (45) is shown. The initial state of Alice and 
Bob is a two mode squeezed state with r = 0.2. Other parameters are A = 300. 7 = 0.002 
(right). The short time behavior is shown in (a). Increase in system-bath coupling and 
acceleration both suppresses the fidelity. The long time behavior in (b) is damped oscillatory 
similarly to the temporal behavior of entanglement. 

evolution of F is plotted. The initial decay (a) is followed by the oscillations (b) with the 
fundamental oscillator frequency fl. We see that the acceleration suppresses fidelity similarly 
to the effect of temperature. 

In the scheme of quantum teleportation for qubits or continuous variables, after Chris's 
state is mixed with Alice's, their states are measured and the measurement result is send 
to Bob in a distance. In the continuous variable quantum teleportation, the homodyne 
measurement on Alice and Chris's state projects their state to the maximally entangled EPR 
state. The imprecision involved in the measurement can be included, for instance, by making 
Gaussian smearing of the EPR state. Here we take into account the finite uncertainty in the 
projected state by considering a two- mode squeezed final state with a squeezing parameter 
r2. We denote the measurement result as a complex number f3 = (3^ + if3i, where f3ji and 
(3i are both real numbers. Thus after the measurement, the original quantum state of Alice 
and Chris is projected to 

00 

\P)ac = cosh-V2Etanh"r2£'(/3)|n)A|n)c, (46) 

n=0 
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where the translation operator D{(3) is given by 



In the hmit r2 ^ oo, \P)ac becomes the maximally-entangled EPR state. 

After the measurement of Alice and Chris's state, the density matrix for Bob becomes 

Pb = Nbac{P\o)cPab c(tt|/3)^c, (48) 

where Nb = cosh(ri — cosh(ri + r2)/7r is a normalization constant. With this normaliza- 
tion, P{l3) = TtePb is equal to the probability to obtain the measurement result f3. Bob's 
state is shifted depending on the measurement result P so that the final state of Bob is 
Pout — D{/3)pbD^{/3). The teleportation fidelity F is an overlap between Bob's final state 
Pout s-nd Chris's initial state \ac) and given by 

B{a\pout\a)B ^^g^ 



Suppose we do not make post-selection following the measurement, the fidelity becomes, 
instead of (49), 

B{a\pB\a)B . ^ 

P{(5) ■ ^ ' 

In this case, the averaged fidelity over all measurement outcome (3 is 

Fav = / (f(5P{(5)F{(5) = j (f(5B{a\pB\a)B. (51) 

Similarly to (46) we write the initial two-mode squeezed state (31) in the Fock space repre- 
sentation as 

oo 

\AB) = cosh-Vi ^tanh'^riln)^!^)^. (52) 

n=0 

In the absence of the environment, from (46), we obtain Fav = (coshri)"^ exp —lap/ cosh^ ri 
in the r2 — oo limit. Then Fav ^ as initial squeezing ri gets larger. This tells us that 
quantum teleportation cannot be successful without post-selection. 

On the other hand, with post selection, the average fidelity can be shown to be 

F„ = l±^. (53) 

For a coherent state ri, Fav = 1/2 and for any positive squeezing, Fav > 1/2. Since Fav = 1/2 
is a maximal value achieved by the classical method[26]. We see that the initial entanglement 
in Alice and Bob's state is essential for quantum teleportation. 
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Now we consider the effect of the environment on the performance of teleportation. In 
our situation, the environment acts continuously on the system so that the final density 
matrix of Bob depends on how we send the measurement result of Alice and Chris to Bob. 
In order to eliminate this ambiguity, here we consider the efi^ect of the environment up to 
the point when Ahce and Chris's state is measured. We introduce the complex variables 
Aj (i = 1, 6) to simplify the expression as follows: 

Ai = {^X)\,-{^X)l^ + 2i{^X)\^, 

A2 = {^X)\, + {^X)l^, 

A3 = {^X)l,-{^X)l, + 2i{^X)l,, 

A4 = {^X)l, + {^X)l,, 

A5 = {/\X)l,-{^X)l, + 2t{{^X)l, + {^X)l,), 

Ae = {^X)l, + {l^X)l, + 2^{{l^X)l,-{^X)l,). (54) 

With a measurement outcome /3, B{(Apout\oi)B and can be calculated straightforwardly. 
Both can be written in the Gaussian form: 

B{oi\pout\oi)B = A^i exp [-V^{a - I3f - f i(a - - 2V2\a - 

and 

P{(3) = Npiex.^[-T^{a- (3f -t^{a-^f -2T4\a- (55) 

The normalization factors and coefficients in exponents are both dependent on the correlation 
matrix and initial condition. As a result, they have the complicated time dependence as 
follows. We will write q = cosh(2rj) and Sj = sinh(2rj) for i = 1, 2 below: 



IN. 



B 



Pf!+pfl + 2qfif2 slA, 



dzdidl d 



2 



r 



p/l/2+p/l/2 + g(|/lP + |/2P) Slc^ 

dsd^di d2 
C3 = C2 + A2, 

C4 = C2 + A2 



61 = AiA^ + AiA^ + 2C3A5A6 



C2-I' 
15 -r ^lAg + 2c3^5/_i6, 

62 = A1A5A6 + AiA5A6 + C3(|A5|2 + lAeH, 
di = (c2 + A2)'-|Ai|2, 

d2 = 4 + C^l(c2 + l)^-2s^(c2 + l)(c2 + A2), 

^3 = q'-\p\', 
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^4 = (di(A4 + l)-62)'-(rflA3-6lf, 

p = A3-d-^[AiA2 + AiA2 + 2c3A5A6], 

- d2-^dr's^[(|Ai|' + C3(c3 + C4))(A^Ai + A^Ai) + 2A5A6(|Ai|2(2c3 + C4) + 0^04)], 

q = A4 + l-dr'[AiA5A6 + AiA5A6 + C3(|A5|2 + |A6|')] + 

dj^dr's^ [ + C3(c3 + C4))(A5AiA6 + AsAiAe) + (|Ai|2(2c3 + C4) + clci){\l^^\^ + lAel")], 

/l = Ml-S^[A5(|Ai|2 + C3C4)+AiA6(c3 + C4)], 

/2 = -s^[A5Ai(c3 + C4) + A6(|Ai 1^ + 0304)], 

r3 = di(0?iA3-6i)/4, 

r4 = cii(di(A4 + l) -62)M, 

^Pi = 47/- (56) 
0(4' 

The averaged fidehty over all measurement outcome (3 is 

(57) 



where 

= Vl-\V^\\ (58) 
In Fig. 10, the time evolution of the averaged fidehty is plotted. For the final state to be in 




Figure 10: The temporal evolution of the fidelity for the two mode squeezed final state with 
the squeezing parameter r2 = 1.0 is shown. The initial state is a two mode squeezed state 
with the squeezing parameter ri = 1.0. 7 = 0.002, and A = 300. The short time behavior in 
(a) shows that strong system-bath couphng and acceleration both suppress the fidelity. The 
long time behavior in (b) is damped oscillatory. 
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a squeezed state with a finite squeezing parameter r2, the initial fidehty is not one. The short 
time decay rate is increased with stronger system-field coupling and with larger acceleration 
of Bob as seen in Fig. 10(a). The long time evolution exhibits damped oscillations. The 
larger the acceleration, the larger the damping rate results. 

The expressions in (56) become much simpler if there is no interaction with the environ- 
ment. In this case, Ai, A3, Ag — > and A2, A4 — > ci, A5 — > —Si. Other terms are 



and 



p,/2,ri,r3,6i 



(59) 



C3 
C4 

d2 

Q 
ds 

di 
fi 

r2 



PI 

di 
b2 

From (55), we have 



and 



Then from (49), 



Cl + C2, 

(Ci + C2)(ci + 1) 

02 -1 

(cl+l)2(c2+l)^ 

2, 
4, 

(C1 + 1)(C2 + 1)' 
(Cl + C2)^ 

(Cl + C2)'(Ci + 1)(C2 + 1) [(Ci + 1)(C2 + 1) + S1S2] , 
1 - 



S1S2 



(C1 + 1)(C2 + 1)' 

Cl + C2 



(C1 + 1)(C2 + 1)' 

4:Nb{Ci+C2) 



(C1 + 1)(C2 + 1)' 
{Ci+C2)\ci + lf{c2 + lf, 
(Ci + C2)sl. 



B{a\pout\o:)B 



B 



(C1 + 1)(C2+1) 



,-2|a-/3pr2 



(C1 + 1)(C2 + 1) 



exp 



-2|q;-/5| 



2 C1C2 



,-2|a-/3|2r4 



S1S2 + 1 



(C1+1)(C2 + 1) 



(60) 



(61) 



(62) 



(63) 
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and 

_ 2 cosh(ri - r2) cosh(ri + , . 

^'"'^ - {c^ + l){c2 + l)-s^S2 ^ 

follow. When the final projected state becomes maximally entangled, r2 oo, then C2, S2 —>■ 
oo and 

= ^ + *";"''•■> . (65) 

Thus we reproduced Eq. (53). Furthermore, when the initial state of Alice and Bob becomes 
maximally entangled, ri oo, then Fav — > 1- Thus there is no loss of fidelity in this limit. 
In the proposed quantum teleportation scheme, the state for Ahce and Chris after 

measurement is assumed to be maximally entangled. It is the Bell state for the qubit 
tclcportation[19] and the EPR state for the continuous variable teleportation[47]. The latter 
can be obtained by taking r2 — >■ oo in (46) as 

oo 

\f3)AC--T.D{f3)\n)A\n)c. (66) 

n=0 

This is the EPR state shifted by the displacement operator D{P). In this limit, Bob's state 
after measurement becomes 

-j^ oo 

Pb^- Yl A{n\c{n\D\(3)\a)cPAB c{a\Dmm)A\m)c. (67) 

^ n,m=0 

The elements in (55) become somewhat simpler as 

B{a\pout\a)B = N2 exp [-T^ia - fif - T^{a - - 2TQ\a - - T^{a - P) - r^{a - P) 
and 

P((3) = Np2exp[-rs{a-(3f-rs(a-P)''-2r^\a-(3\ 
The coefficients are 



No = 



d, = (A4+i)2-iA3r, 

J 2 I |2 

de = r - \s\ , 

r = (A2 + l)rf5-(A4 + l)(|A5p + |A6|2), 

s = Airfs -2(A4 + 1)A5A6, 

u = 4 - (A4 + 1)A5, 

V = -(A4 + 1)A6, 
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Tio 



1 



44 
1 

44 
1 



sii^ + st;^ + 2ruv 



SUV + SUV + r{\u\'^ + \v\'^) 



+ 



A4 + I 



— [suw + svw + ruw + rvw\ + —, 

44 4 
= A34-\ 

= (A4 + l)4-\ 

1 



444 



sw + sw +2r\w\ 



(68) 



The averaged fidelity over all measurement outcome (3 in this case is 



d'^/3B{a\pout\a)B 



d'pPi(3)Fi(3) 



2d^ 



1/2 



(69) 



where 



Til 



q-|r5l 



= 4^(r5r? + r5r2-2r6|r7|2). 



(70) 



In Fig. 11, the time evolution of the average fidelity for the maximally entangled final state 
is plotted. For this final state, the initial average fidelity reaches the unit value as initial 
squeezing becomes infinitely large as seen in Fig. 11(a). The long time evolution is again 
damped oscillatory with the damping rate enhanced by acceleration as shown in Fig. 11(b). 
In the absence of the environment. 



a,q,r, Fg, Fy, Fg, Fio, Fn 

Other coefficients will be simplified accordingly as 

2 

n{c^ + 1) ' 

-51 



0. 



(71) 



4 
4 



Ci + l' 

(ci + 1)^ 
4(ci + 1)^ 
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Figure 11: The temporal evolution of the average fidelity for the maximally entangled final 
state is shown. The initial condition is a two mode squeezed state with the squeezing 
parameter ri. The short time regime is shown in (a) with a — 0. The long time regime is 
shown in (b) with ri = 1.0. Other parameters are 7 = 0.002, A = 300. 



Thus we have 



and 



r 
u 

NP2 



Si 



ci + L 

2(ci + l)^ 

2(ci + l)(ci + si + l), 
2 

7r(ci + l)' 
1 



Ci + 1 



B{a\pout\a)B 



1 



7r(ci + 1) 



The fidelity becomes 



F — exp 



-2\a-(5\^\l 



Sl 



Ci + 1 



(72) 



(73) 



(74) 



(75) 



(73), (74), (75) can also be directly obtained from (63) by taking the limit of r2 — * 00. The 
same expression of Fav as in (65) follows, consistent with the previous result for a generic 
initial and final squeezed state. In particular, in the absence of the environment, we achieve 
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the unit fidelity for the quantum teleportation with the initial and final EPR state, which 
is the case originally studied in [47]. We now consider the effect of the environment in this 
special case. Bob's state after measurement can be written similarly as in (48) as 



Pb = - Y] A{n\c{n\D\P)pABcD{(3)\m)A\m)c- 



(76) 



n,m=0 



Here we consider the influence of environment on all three observer's states and wrote their 
state in the density matrix form as Pabc- We obtain the averaged fidelity as 

1 



TP — 
^ av — 



;Sn + l)(E22 + l) 



E2 



12 



1/2 ' 



(77) 



where are components of a 2 x 2 matrix made up of the components of the 6x6 matrix 
S in (21). We first write S as 
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Then 



(79) 



In Fig. 12, the time evolution of the average fidelity for the maximally entangled initial and 
final state is plotted. The decay is faster for larger acceleration. After the damped oscillatory 
behavior, all curves appear to reach nonvanishing asymptotic values. We will now show this 
is indeed the case. 

In the long time limit, the state of Bob in the absence of Alice and Chris will thermalize 
at the Unruh temperature Tu similarly to the single particle quantum brownian motion at 
finite temperatures [48]. On the other hand, the state of Alice and Chris in the absence of 
Bob will be in the steady oscillatory state (see the discussion in [25]). Now let us first assume 
the interaction of Bob with the environment is negligible. In this case, E^ = 0. In the weak 
coupling limit [2 5], 




Then 




(80) 



(81) 
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Figure 12: The temporal evolution of the fidelity for the maximally entangled initial and 
final state is shown. 7 = 0.001, A = 50. 



follows. This will give the averaged fidehty 



- 0.70. 

V2 



(82) 



Thus although there is still information loss, a part of quantum information manages to 
escape making quantum tclcportation successful [46]. This case can be compared to [15], 
where the authors pointed out that the unitary quantum gates between the original qubit 
with one of the Bell pair prior to the measurement ruin the unitarity of the transfer matrix, 
and to [16], where the author included random unitary gates in quantum circuits estimated 
the fidelity to be 0.85. 

Now we take into account the interaction between Bob and the environment. This will 
add the thermal factor in the induced fluctuations of Bob's state in the long time limit as 



coth(Q/27V) ( 1 
4 I 1 



Then 
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The averaged fidelity in this case will decrease to 

2 



(coth((]/2Tt;) + 4)V2(coth(l]/2T[;) + 8)1/2 



< 



45 



0.30 < 0.5, 



(83) 



(84) 



(85) 



where the equality holds for zero acceleration {Tjj = 0). We see that quantum teleportation 
fails in this case. In this sense, quantum information will be totally lost. 
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5 Conclusion 



We studied the quantum entanglement dynamics of static and accelerated observers. Due 
to the interaction with the environment, quantum states of the observer suffer decoherence 
causing their entanglement to decay. In addition, due to accelerated motion, the observer 
sees the vacuum as a thermal state, which enhances the rate of disentanglement in observers' 
quantum states moving in the vacuum. 

We derived uncertainties for the quantum states of two and three observers and saw 
their dynamical evolution. Since our calculations are exact, uncertainty relations are always 
satisfied as they should be. They are no longer guaranteed under perturbative approxi- 
mations. The evolution of uncertainties shows damped oscillations whose amplitudes are 
enhanced by the Unruh effect. The larger the acceleration of the observer, the larger the 
oscillation of uncertainties because the effective temperature (Unruh temperature) for the 
thermal fluctuations is proportional to the acceleration. 

The entanglement dynamics between two observers, one static and the other under con- 
stant acceleration was analyzed. From the uncertainty of the partially transposed density 
matrix, the negativity and the log- negativity were calculated. They showed that the Unruh 
effect on Bob dynamically enhances the uncertainty after the partial transpose leading to 
the suppression of bipartite entanglement between Alice and Bob' state. 

We further studied tripartite entanglement between two static and one accelerated ob- 
server. The acceleration modifies the three-mode uncertainty similarly to the two-mode 
uncertainty. The three-mode uncertainty shows the damped oscillations while satisfying the 
uncertainty relation. The environment dynamically modifies initial entanglement between 
Alice and Bob and induces entanglement between Alice and Chris. As a result, bipartite en- 
tanglement Eab and Eac change in time in an opposite manner. No bipartite entanglement 
is induced between Bob and Chris. Bipartite entanglement is suppressed by acceleration of 
observers. From the combinations of bipartite entanglement among three parties, the gen- 
uine tripartite entanglement Eabc was calculated. The positiveness of Eabc indicates that 
the monogamy inequality is always satisfied. For our initial condition, a product state of a 
two-mode squeezed state of Alice and Bob and a coherent state of Chris, Eabc vanishes ini- 
tially. It becomes positive due to the effective interaction between Alice and Chris induced 
from the vaciium at intermediate time scale. In a meantime, decoherence induced from 
the environment causes Eabc to decay at long times. Similarly to the bipartite entangle- 
ment dynamics[25], the environment plays a dual role. It induces multipartite entanglement 
among subsystems but it also is the source of decoherence and dissipation that decrease 
entanglement causing information loss. 

Next we studied quantum teleportation protocols between two static and one accelerated 
observer. We first looked at the case with arbitrary final states of Alice and Bob by summing 
over all possible final states of them. What we obtained is an one-mode Q function for Bob's 
state. In the long time and weak coupling limit for a zero acceleration, we saw that the 
fidehty reaches one and there is no information loss in this case. Then we considered the 
generalized final state to which the quantum state of Alice and Chris projected after the 
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measurement to be an arbitrary two- mode squeezed state with the squeezing parameter r2. 
The timc-dcpcndcnt tclcportation fidehty is calculated for the measure of success of quantum 
teleportation. We showed that after performing the proper unitary transformation on Bob's 
quantum state based on the measurement result sent to Bob, the fidelity is recovered to a 
unit value in the limit of ri,r2 oo in the absence of the environment. For the generic 
value of ri,r2, however, the fidelity is always smaller than 1. Thus there is information loss 
even in the absence of the interaction. Furthermore, even if the initial state of Alice and 
Bob is maximally entangled, the fidelity is still smaller than 1 for generic finite values of r2. 

In the presence of the interaction with the vacuum, the fidelity shows damped oscilla- 
tions similarly to entanglement. This indicates the close relation between entanglement and 
fidelity. In the absence of environment, there is a direct relation between fidelity and entan- 
glement for maximally entangled final states as seen in (53). With the generic final state 
in the presence of the environment, the relation is not so straightforward. The damping of 
fidelity is enhanced by the acceleration. The Unruh effect overall suppresses fidelity causing 
more information loss. We obtained the fidelity expression for the maximally entangled final 
state case r2 — > oo in (69). This type of final states is usually assumed in current experiments 
for quantum teleportation [49]. We saw that the temporal behavior of the fidelity consists of 
damped oscillations similarly to the case with generic r2. The maximum value of the fidelity 
is achieved only at the initial time, which becomes a unit value as ri goes to infinity. 

We also examined the case when the initial and the final state are both maximally 
entangled (ri,r2 — > oo). In this case, the teleportation is perfect (the fidelity reaches the 
unit value) in the absence of the environment or other sources of information loss. In the 
presence of the interaction with vacuum, however, the fidelity is shown to decay. Taking the 
long time limit, we obtained the asymptotic value of the fidelity. If we ignore the interaction 
of Bob's state with the vacuum, we obtain F^v ~ 0.7 and quantum teleportation can still 
be viewed as successful. Proper encoding of the message with quantum error correction 
may help retrieving the complete information. If we take into account the interaction with 
vacuum for all three parties, we obtain Fav ~ 0.3, thus quantum teleportation is no longer 
successful. 

Implications of our results to black hole information problem can be summarized as 
follows. There are two fundamental sources of nonunitary evolutions and information loss. 
(1) the boundary condition of the evolution of the total system (intrinsic loss). (2) the 
interaction of the relevant quantum mechanical system with the surrounding quantum field 
(induced loss). (1) can be attributed in the laboratory to the special choice of boundary 
conditions imposed by measurement. (2) is present in any spacetime and responsible for the 
quantum/classical correspondence of quantum systems of large degrees of freedom [43]. In a 
black hole spacetime, the singularity provides an effective boundary on the quantum state of 
infalling matter and the incoming Hawking particle. The origin of nonunitary evolution in 
(1) may be explained by the remnants or wormholes[12]. (2) is also present as the interaction 
among infalling matter. Hawking particles, and gravitational field of a black hole. Thus in 
addition to the correlation in a Hawking pair, there is the induced three-body correlation and 
tripartite entanglement between the matter and a Hawking pair, which vary dynamically due 
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to gravitational interaction. In the practical setting of the laboratory, there are many other 
phenomenological loss of information such as photon loss during the optical transmission, 
leakage, etc. We did not consider these effects in this work. 

At the present stage without complete knowledge of quantum gravity, we do not know 
details about the final state of the black hole. Assuming a maximally entangled final state, 
we may retrieve quantum information from a black hole just as in quantum teleportation as 
suggested in [14] provided that (I) the induced loss due to vacuum is negligible. Since the 
environment is always present, this type of loss is unavoidable. (11) The proper postselection 
is performed. This can be viewed as an extension of the time-neutral formulation of quantum 
mechanics[50] applied to spacetime physics, in which both initial and final states are specified. 
Decoherence history interpretation of the formulation in view of cosmology is given in [51]. 
In the meantime, performing the postselection in the laboratory is now very common in 
quantum information science. We hope our attempts using accelerated observers provide a 
practical setting to understand the black hole information problem. 
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